SCALING EXPONENT FOR THE HOPF-COLE SOLUTION 
OF KPZ/STOCHASTIC BURGERS 



M. BALAZS, J. QUASTEL, AND T. SEPPALAINEN 

Abstract. We consider the stochastic heat equation 

dtZ = d^Z-ZW 

on the real line, where W is space-time white noise. h{t,x) — — log Z{t,x) 
is interpreted as a solution of the KPZ equation, and u{t, x) — dxh(t, x) as a 
solution of the stochastic Burgers equation. We take Z(0, x) = exp{i?(a;)} 
where B{x) is a two-sided Brownian motion, corresponding to the station- 
ary solution of the stochastic Burgers equation. We show that there exist 
< Ci < C2 < oo such that 

cif/^ < Var(logZ(t,a;)) < C2t^^^. 

Analogous results are obtained for some moments of the correlation func- 
tions of u(t, x). In particular, it is shown that the excess diffusivity satisfies 

cii^/' < D{t) < C2t^/\ 

The proof uses approximation by weakly asymmetric simple exclusion pro- 
cesses, for which we obtain the microscopic analogies of the results by 
coupling. 



1. Introduction 

1.1. Physical background. The Kardar-Parisi- Zhang (KPZ) equation [13] 
is a formal stochastic partial differential equation for a random function h{t, x), 
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t>0, X eR, 

dth = -X{d^hf + udlh + aW (1.1) 

where z/ > and cr, A 7^ are fixed parameters and W{t,x) is Gaussian space- 
time white noise 

E[Wit,x)W{s,y)] = 6{t-s)6{y~x). 

It is widely studied in physics as a model of randomly growing interfaces. The 
derivative u = d^h should satisfy the stochastic Burgers equation, 

dtu = -Xd^u'^ + udlu + ad^W. (1.2) 

Using renormalization group methods physicists have computed the dy- 
namic scaling exponent ([ID], [13], [1]) 

z = 3/2. 

Roughly, this means that one expects non-trivial behavior under the rescaling 

h,{t,x) =e'/^h{e-''/\e-^x). 

For the totally asymmetric exclusion process and the polynuclear growth mod- 
els, which can be thought of as discretizations of (11. ip . it is now known rigor- 
ously [9j, |20j that in a weak sense, 

Var(/i,(t, x)) ~ t^/'gscit-'^Hx - vt)) (1-3) 

for an explicit v and scaling function Qsc related to the Tracy- Widom distri- 
bution. Note that these models are in some sense exactly solvable. 

(11.11) and (11.21) are ill-posed because the quadratic non-linear term cannot 
possibly make sense for a typical realization, which, in the case of (II. ip is 
expected to look, in x, locally, like a Brownian motion with variance z/~^cr^. 
Formally applying the Hopf-Cole transformation 

Z{t, x) = exp{-Xi^-^h{t, x)} (1.4) 

to (II. ip leads to the stochastic heat equation 

dtZ = vdlZ - Xiy-^aZW. (1.5) 

The advantage is that (II. 5p is well-posed [23]. We do not attempt to justify 
the manipulations leading to (11.51) . We define h and u = d^h through (II. 4p . 
These Hopf-Cole solutions are expected to be the physically relevant solutions 
of (dH) and (ra . 

Z can be thought of as an asymptotic model of a directed polymer. There 
is a Feynman-Kac formula 

Z{t,x) = E^[:exp{-(3 [ W{s,b{s))ds}: Z{0,b{t))] (1.6) 

Jo 
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where the expectation is over an independent Brownian motion b{s), s > 
starting at x, of variance u, j3 = Xu^^a, and :exp: is the Wick-ordered 
exponential (see [18] for details). The reason we write (11. 6p is to draw attention 
to the analogy with directed polymers, a typical model being 

z{n,x) = E4exp{-l3Y.'L=iXim,Sm)}] (1-7) 

where X(m, r), m G {1,2,...}, r G {. . . , — 1, 0, 1, . . .} are independent and 
identically distributed random variables, and Sm is a simple random walk 
starting at x. Assuming reasonable decay on the tails of the X's, it is expected 
[14j that for any /?, 

Var(log2(n,0)) ~ cn^^ with x = 1/3- (1-8) 

Little is known rigorously. Bounds analogous to x ^ [3/10, 1/2] are obtained 
in [0], [H], [IS], [IH], and [22] • The closest results with x = 1/3 are those of 
[H [12] for certain last passage percolation models, which are obtained in the 
/9 — *• oo (zero-temperature) limit. Note the contrast with dimensions d > 3 
where the polymer is known to be diffusive for small /3. 

1.2. Mathematical background. We now survey what is known rigorously 
about (II. ip . In terms of well-posedness, the technology at the present time [21j 
can only handle far smoother noise terms than the white noise. An unusual 
type of Wick product version of the problem has been introduced [TT]. But 
besides requiring fairly smooth noises, this does not have the scaling expected 
[8], and is therefore believed not to be physically relevant. 

The idea in [5], which leads to what appears to be the physically relevant 
solution, is to smooth out the white noise in space a little, and then use the 
Hopf-Cole transformation and the tractability of (11.50 to remove the cutoff. As 
this is done, one finds one has to subtract a large constant from the equation. 
The resulting Hopf-Cole solution of (11.11) is given explicitly as the logarithm 
of the well-defined solution of (II. 5p . We now recall the details. 

Let W(t), t > 0, be the cyhndrical Wiener process, i.e. the continuous 
Gaussian process taking values in Hy^J (M) = na<-i/2H^^{R) with 

E[{^, W{t)){i;, Wis))] = min(t, s){^, i^) (1.9) 

for any (fjip & C^(M), the smooth functions with compact support in M. Here 
H^^{R), a < 0, consists of distributions / such that for any ip G C^(]R), ipf 
is in the standard Sobolev space if~"(]R), i.e. the dual of if"(M) under the 
pairing. iJ~"(R) is the closure of C^(M) under the norm J (H-|t|~^")|/(t)p(it 
where / denotes the Fourier transform. 

The distributional time derivative W{t,x) is space-time white noise. These 
Sobolev spaces are the natural home, as can be seen by restricting to a finite 
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box and writing a Fourier series for white noise, with independent and iden- 
tically distributed Gaussian variables as Fourier coefficients. Note the mild 
abuse of notation for the sake of clarity, as we write W{t, x) even though it 
is a distribution on [t, x) G [0, oo) x M as opposed to a classical function of t 
and X. 

Let J-'{t), t > 0, be the natural ffitration, i.e. the smallest a-field with 
respect to which W{s) are measurable for all < s < t. Let Q G C^(M) be 
an even, non- negative function with total integral 1 and for k, > 0, Q^{x) = 
K~^Q{n~^x). The mollified Wiener process is given by x) = {T-xGK,W(t)) 

where T-xf{y) = f{y + x) is the shift operator. The distributional time deriv- 
ative of W^ltjx) is the Gaussian field with covariance 

C,{x - y)6{t -s)= ElW'^it, xW'^is, y)] (1.10) 

where 

= / g{^x - y)g{-'^x - y)dy, CM = ^''WOWl 

Assume the initial data h{0, x) is a random continuous function on M such 
that for each p > there is an a = a{p) < oo such that 

sup^^^e-''\''\E[e-P''^''^] < oo. (1.11) 

The mollified KPZ equation is 

dth"" = -A[(9^/i^)2 - C«(0)] + udlh^ + aW\ (1.12) 

The following summarizes previous results, mostly from [5]. 

Proposition 1.1. Let h{0,x) be a random initial continuous function satis- 
fying lil.ll]) and independent of the white noise W. 

1. For each /t > 0, there exists a unique continuous Markov process hf^{t), with 
probability distribution on C([0, T], C(M)), adapted to J-'{t) = a{hQ, Wg, s < 
t), t > 0, and solving U.l^) . 

2. The process Z'^{t,x) = expl—Xu^^h^^tjx)} is the unique adapted (mild) 
solution of the ltd equation, 

dtZ" = udlZ" + Az/" VZ'^M/^ Z'^iO, x) = exp{-Az/"^/i(0, x)} (1.13) 

and has the representation 

Z^{t, x) = E4exp{-f3 f W^{s, b{s))ds - ^a(0)t}Z'^(0, bit))]. 

Jo ^ 

3. Z'^(t,x) — >• Z(t,x) almost surely, uniformly on compact sets of [0, oo) x 
M where Z is the unique adapted (mild) solution of U.5\) with Z{0,x) = 
exp{—Xh'~^h{0,x)}. Furthermore, Z(t,x) > for all t > 0, x & M. almost 
surely [T7] . 
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4- The Pk are tight as measures on C{[0,T],C(R)). The limit process h{t) 
coincides with 

h{t, x) = -A" V log x). (1.14) 

5. Define 

u = d^h (1.15) 

in the sense of distributions. If we start Z with initial data Z{0, x) = exp{B{x)} 
where B{x) is a two sided Brownian motion independent of W with vari- 
ance z/~^o"^, then u is stationary in both space and time. In this sense, 
Gaussian white noise with variance v~^a'^ is invariant for U.^) . Station- 
arity here means stationarity of {Tx(f,u{t)) for smooth functions of compact 
support ip{x) where {Trc{p){y) = {p{y — x). The corresponding h(t,x) and 
Z{t,x) are not stationary in time with these initial data, but the increments 
Dsh{t, x) = h{t, X + 5) — h{t, x) are space and time stationary. 

6. Let h be the Hopf-Cole solution of M.l\] . as in Ij^l.l4 ). Then 



h^{t, x) = -f'^hij-H, j-^x) (1.16) 

is the Hopf-Cole solution of U.l\) with new coefficients 

A^ = ^-/3+2^, j,^ = ^-P+2^^ a^ = ^-^^a. (1.17) 

If we start 1^1. S\} in equilibrium, ie. u{0) is a white noise with variance v^^a^, 
then the time reversed process u(T — t),tE [0,T) is a solution of U.^) with 

A replaced by —X, and the spatially reversed process h{t, —x) h{t,x). 

Note that [5J only consider the case X = u = 1/2 and a = 1, but their 
proofs work in general. 6 is not stated there, but it follows readily from their 
methods. 

The Hopf-Cole solution fll.l4p of KPZ (11. II) is obtained as a limit of solutions 
of (I1.12p . i.e. after subtraction of a divergent term (7^(0). An important open 
problem is to show that a corresponding version of (11.11) with an appropriate 
Wick ordered nonlinearity is well-posed. We do not address this issue here. 
Since (11.141) is expected to be the relevant solution, we simply study it directly. 

1.3. Statement of results. We can now state our main results about h. 

Theorem 1.2. Let h{t,x) be the Hopf-Cole solution of U.l\) as in (TJ^ 
with Z{t,x) the solution of U.5\) . with initial data Z{0,x) = exp{i?(a;)} where 
B{x) is a two sided Brownian motion independent ofW with variance v~^a'^. 
Let Var(/i(t, x)) denote the variance of h{t,x). Var(/i(t, x)) is a symmetric 
function of x, non- decreasing in \x\, and 

Var(/i(t,a;)) - |x| > 0. (1.18) 
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Furthermore, there exist cq = Co(cr, z/, A) < oo, and Ci = Ci{a,h',X) < oo, 
C2 = C2(m, 0", i^, A) < 00, such that for t > cq we have 

Cxt^l^ < Var(/i(t, 0)) < C^t^/^ (1.19) 

and, for 1 < m < 3, 

< f |^|m-2 [Var(/i(t, x)) - dx < Cat^'"/^ (1.20) 



Remark 1.3. The dependence of the constants Cq and C on m, a, z^, A is as 
follows: We can take 



Co = coa-'u'X-\ C, = Ciaf^+V-f-^A^, = C{m)a'^+'u-f~'X'f 

for some co,Ci G (0, 00) and C{m) = ^(3_^(^_i) , 1 < m < 3 and C(l) = 
C/4. Here m = 1 refers to (11.191) . The lower bound holds for all m > 1. 

Remark 1.4. Theorem 11.21 tells us that the scaling exponent for the quantities 
studied here follows the physical prediction, as can be seen by integrating 
(11.31) . This provides considerable support for the notion that this process h is 
the sought after solution of (11.11) . even though is not known presently how to 
show directly that h{t,x) = — X'^u log Z{t,x) solves (II. ip . or indeed, what it 
means to solve (II. ip . 

Remark 1.5. Sometimes we want to indicate the dependence of the solution on 
the parameters by writing h(t, x; A, u, \fv). It is interesting to take A = 1, = 
V because hit^ x; 1, u, ^/u) = vh{v~^t, 1, 1, 1) corresponds to the solution 
of dth = —{dxhY + z/9^/i + ^/vW and note that the result implies that there is 
a z/q > and fixed < ci < C2 < 00 such that cit^/^ < Yai{h{t, 0; 1, z/, y/v)) < 
C2t^^^ for t > 1 for all v <vq. 

Now we turn to results about the correlations of u(t,x). Throughout we 
will assume that it is in equilibrium with initial data white noise with variance 
z/~^(j^. By definition {(p,u{t)) = — J (p'{x)h{t, x)dx for (p G C^(M). For each 
fixed t > 0, define a bilinear functional on C^(M) x C^(M) 

Bt{^,^) = E[{^,u{t)){^,u{0))]. 

The following proposition provides us with our definition of the space-time 
correlation measure of the stochastic Burgers equation. 

Proposition 1.6. For each t > there is a unique probability measure S{t, dx) 
on M such that for (p,ip E 



S{t,dx). (1.21) 
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S(t, ■) is symmetric: S(t,A) = S(t, —A) for Borel sets A where —A = {—x : 
X G A}. The connection with the process h is that 



Var(/i(t,x)) = x + 2 / {z - x) S{t,dz) (1.22) 

Jx 

for a; G M and t > 0. 

The route to constructing S(t) and proving (11.221) is somewhat circuitous. 
The measure S{t) is constructed as a weak hmit from the rescaled correlations 
of a particle process. Then we show that, in the sense of distributions, 

Idl VaT{h{t, x)) = S{t, dx). (1.23) 

Finally, after studying solutions of the stochastic heat equation, we can deduce 
(ll.22p . Here are the bounds on S{t). 

Theorem 1.7. Let u{t) = dxh{t) be the distributional derivative of the Hopf- 
Cole solution fll.l4p of (11.11) . in equilibrium with initial data white noise with 
variance z/~^cr^ and let S{t, dx) be the space-time correlation measure defined 
through (ll.2ip . With the same constants as in Theorem for t > cq we 
have 

^^^2m/3 < j < c-st^Ws^ 24) 

The upper bound holds for 1 < m < 3. The lower bound holds for all m > 1. 
In particular, the diffusivity 

1 /■„2, 



D{t) = - J x'S{t,dx) (1.25) 
satisfies 

Cit^''^ < D{t) < C2t^'^. (1.26) 

We turn to proofs. The first issue is to develop the connection with the 
exclusion process. 

2. Weakly asymmetric simple exclusion 

We consider nearest neighbour (i.e. simple) exclusion on Z (ASEP) with 
particles attempting jumps to the right at rate p = 1/2 and to the left at rate 
g = 1/2 + e^/^ with e G (0, 1/4). This is a system of continuous time random 
walks jumping to the right at rate p and to the left at rate g, with the rule 
that jumps to already occupied sites are not realized. Hence the occupation 
variable can be taken to be ri{t, x) = 1 or depending on whether or not there 
is a particle at x & 1^ at time t. One of the most important properties of this 
system is that it preserves Bernoulli product measures with any density p G 
[0, 1]. Here we take p = 1/2, i.e. we take as initial configuration independent 
Bernoulli {r]{0,x)}, x G Z, with density 1/2. 
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Let t) = 27] — 1 and define the height function 

Ut.x) = \ -2iV(i,0), x = 0, (2.1) 

I -E.<,<o^(^'2/)-2iV(^'0)' ^<0' 

where N{t, 0) is the current across the bond (0, 1) up to time t, i.e. the number 
of particles that jump from to 1 minus the number of particles that jump 
from 1 to in the time interval [0, t]. For a; G M and t > let h^{t, x) denote 
the rescaled height function; 

K{t, x) = e^/^ [U^'h, [e-\]) - Vet) (2.2) 

where = \^~'^^'^ — ""^^^ and the closest integer \x\ is given by 

[x] = Lx + |J. (2.3) 

We think of as an element of the space -D([0, oo); i5^j(M)) where D refers 
to right continuous paths with left limits. -Du(M) indicates that in space these 
functions are equipped with the topology of uniform convergence on compact 
sets. Because the discontinuities of heit, ■) are restricted to e(l/2 + Z), it is 
measurable as a D„(M)-valued random function (see Sec. 18 of [7j.) Since the 
jumps of h^it^ ■) are uniformly small, local uniform convergence works for us 
just as well the standard Shorohod topology. The probability distribution of 
the process on -D([0, oo); D„(M)) will be denoted V^- 

Proposition 2.1. [5J yls e \ 0, the distributions Vg converge weakly to V, the 
distribution of the Hopf-Cole solution h f ll.Mp of (11.11) vuith X = 1/2, v = 1/2, 
and (7 = 1. 

Proof. This was proved in [5] using the slightly different height function C^'~'(^5 x) 
related to ours by 

C^«(t, a:) =Ce(t,x)- 2(1 -7^(0,0)). (2.4) 

The result follows for (i;(t,x) because the difference is bounded. Note also 
that [5] makes the height functions continuous in space by linear interpolation 
and uses the smaller path space D([0, cxo); C(M)). This makes no difference 
because h^'^{t, x) and its continuous version are uniformly within distance e^^^ 
of each other. I 

The rescaled velocity field is 

Ue{t, x) = 6-^/^fi{6-h, [e'^x]) (2.5) 

with [■] defined by (12. 3p . The rescaled space-time correlation functions are 
given by 

S,{t,x) = E,[u,{t,x)u,{0,0)]. (2.6) 



For X G 
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let US define a discrete Laplacian by 

A,/(a;) = \e-\i\x + e) - 2f{x) + f{x - e)) (2.7) 
and a discrete absolute value in terms of the closest integer function by 



\x\ 



e e X 



(2.8) 



We begin by building on some well-known properties. The assumption that 
density is 1/2 is used repeatedly. 

Proposition 2.2. 1. For a fixed t > 0, S'e(t,x) is a probability density on M 
and symmetric in x except at x & e{l/2 + Z) . 

2. Seit,x) = AeVaY{hs(t,x)). YeiT{hs(t,x)) — \x\e is symmetric in x ^ 
e{l/2 + Z), non- decreasing in \x\e, nonnegative, and for each fixed {e,t) has 
exponentially decaying tails in x. 

3. For 1 < m < oo, 



x\'^Se{t,x)dx = / A,{\x\'^)[VaT{h,{t,x))-\x\,]dx. 



(2.9) 



Proof. To see that Si;{t,x) is a probability density, use the well-known con- 
nection with the second class particle: 



S,{t,x)=e-'Pl/'Me~H) = [e''x]}. 



(2.10) 



1/2 

Here Pe is the coupling measure of two ASEP's that start with one discrep- 
ancy at the origin and Bernoulli (1/2) occupations elsewhere, and x(-) is the 
position of the second class particle. This setting is discussed in Section HI A 
proof of fl2.10p can be found for example in Symmetry of S^{t,x) can be 
seen from the definition 02.51) and the fact that 

?j{t,k) = l-r]{t,-k) (2.11) 



defines an ASEP rj equal in distribution to rj. See [20] for the explicit compu- 
tation that proves 2. 

We now work towards 3. Let N{t, x) denote the current across the bond 
between site x and x + 1 up to time t. We start by checking that 



Gov 



0. 



(2.12) 



y=—x+l 

With rj as in (12.111) . an ?7-particle jump from x to y is the same as an ?7-hole 
jump from —x to —y. Hence N{t, 0) = N{t, —1). By the distributional equal- 



ity v = V, Gov N{t,0),Z:=-x+i^it^y) 



Gov 



Gov 



is verified. 



Gov 



N{t,0),E:=-x+Mt,y) andim 
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Combining Lemma 3 of [22j with fl2.12p gives 

Vai{Ceit,x))-\x\=ACoY{N{t,0),N{t,x)), x E Z. (2.13) 

The right-hand side of fl2.13p is symmetric in x by invariance under spatial 
translations. 

Next, note that by the finite range of ASEP, for fixed e > and t > there 
exist Ci < oo, C2 > such that 

\CoY{N{t,0),N{t,x))\ < Ciexp{-C2|a;|}. (2.14) 

(see Lemma 4 of [22] for details.) 

We next argue the nonnegativity of fl2.13p . By symmetry it suffices to 
consider x G For x > 0, v{x) = Var((^£(t, x)) and v{x) — have the same 
discrete Laplacian. By (12.131) and (12.141) both v{x) — \x\ and Aiv{x) decay 
exponentially, and, by 1 and 2, Aiv{x) > 0. Then for x G 



00 00 

v(x) — \x\ 



J]^2Aif(£ + 1) > 0. (2.15) 

k=x l=k 

This also shows that v{x) — \x\ is strictly decreasing for a; G Z_|_, and thereby 
(I2.13P is strictly decreasing in 
Putting the scaling into (12.131) gives 

Yax{he{t,x)) - \x\e = AeCoY{N{e-h,Q),N{e-h,[e-^x])). (2.16) 

This proves (3). To prove (4), start with the observation 

-N 



\x\^Ae{\x\e)dx = Q. 

■N 



Then by (2) and by integration by parts (that is, by shifting the integration 
variable) , 

/N f-N 
|x|™^e(t, x)dx = / A,[Vai{h,{t, x)) - \x\,] dx 

■N J-N 



N 



A,(|a;|r) [Var(/i,(t, x)) - \x\e] dx + B.^tAN) 



-N 



where Be^t,m{N) are sums of integrals of \x ± £|^[Var(/i£(t, a;)) — \x\e] over 
intervals of length e around ±A^. By 3 these are exponentially small in as 
£, t and m are fixed. Taking N ^ 00 gives 4- i 

The key technical estimate which will be proved in Section H] is 

Theorem 2.3. With the same constants as in Theorem M.^ for all < e < 
1/4, 1 < m < 3, and t > cq, 

Cit2-/3< j \x\^Se{t,x)dx <C2t^'^/\ (2.17) 
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Corollary 2.4. 1. For < £ < 1/4 and t > cq, 

Cit^/s < Var(/i,(t, 0)) < C2t^''\ (2.18) 

2. For each t > 0, the family of probability measures {5'e(t, x)(ix}o<e<i/4 is 
tight. 

Proof. Part 1 follows from Theorem 12.31 and case m = 1 of of Proposition 
12.21 because A^dxy = for x G [— £:/2,£:/2) and vanishes elsewhere. 

For t > Cq tightness of {5'e(t, x)(ix}o<e<i/4 follows from the upper bound 
in (12.171) . For < t < Cq recall the second class particle connection (I2.10p . 
Proposition 4 in ^2] proves that the second moment 

j \x\lSe{t,x)dx = e^Y.{\-^{e-H)\^) 

is monotone nondecreasing in t. Thus the large-t bound gives the tightness 
for all t > 0. I 



3. Proofs of the main results 

As a preliminary point we discuss the regularity of Var(/i(t, x)). The control 
comes from the weak limit /ig — > h. We have 

1/2 

Var(/i^(t,0)) 



As mentioned in the proof of Corollary 12.41 this last quantity is nondecreasing 
in t. Consequently by the upper bound in fl2.17p and the i.i.d. mean zero 
spatial increments of h^it.x) [see (12. ip ] we conclude that Var(/;,e(x, t)) is lo- 
cally bounded as a function of {t,x), uniformly in e > 0. By the weak limit 
Var(/i(t, x)) < lim^^o Var(/i£(t, x)) and so Var(/i(x, t)) is locally bounded. The 
x-symmetry of Var(/i(t, x)) follows from part 6 of Proposition ll.il or from the 
weak limit —>■ h and the distributional symmetry of Celt, ■). For any fixed 
Xq, h{t,x) — h{t,Xo) is a Brownian motion in x and hence the continuity of 
X I— > Var(/i(t, x)). By studying the stochastic heat equation we prove in the 
Appendix that 

Var(/i(t, x)) — |a;| — » as |x| oo. (3.1) 
We turn to proving the main results. 

Proof of M.21\) . From the definitions we have for test functions v^, G C^(M) 
E[{^\ K{t)){ij', HM)] = E[{¥>, Ue{t)){tp, uM)] + 0{e) 



/y + x\ , (y — x\ ] 
'^[-^)i'y-^)dy\Se{t,x) dx + 0{e). 
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Let S{t, dx) denote a weak limit point of ^^(t, x)dx as £ \ 0. Taking the limit 
in fl3.2l) . the last expression becomes the right-hand side of fll.211) . 

Convergence of the first expectation in (13.21) to the left-hand side of (11.2ip 
follows from the weak convergence ^ h with the following additional justifi- 
cation. Since the limit process h is continuous (in time), the pair (/ie(0), h^{t)) 
converges weakly to {h{0),h{t)). Integration against a C^(M) function is a 
continuous function on Du{M.). By an application of uniform integrability and 
Schwarz inequality, for the claimed convergence it is enough to show the 
boundedness of \ {<f', /^^(t))!^. To see this we first transform the inner product 
(this is the beginning of the computation that one performs to check (13. 2p ): 



= -e'/'Y.^{{k-\)e)1i{e~\k). 



(The constant term on the right-hand side of definition (12.21) vanishes since 
we are integrating the height function against a derivative.) This is a sum of 
independent mean zero random variables, and 

E\{^\ KXm' < Ce' Y.^{{k- \)eY + Ce' J2^{{k- 1)^) V((^ - 

k k,e 

which is bounded uniformly in £. I 

Proof of Theorem |i. 7| and the upper bound of Theorem li.M Note first of all that 



from part 6 of Proposition II. it suffices to prove all results with A = 1/2, 
u = 1/2, a = 1. 

The upper bounds of (11.191) and (11.241) follow from the weak convergence 
and from the upper bounds in (I2.17P and in 1 of Corollary 12. 4[ 

Let 1 < m < 3. For the upper bound of (ll.20p we collect these ingredients: 
Inequality Var(^£(t, x)) — \x\s > from 3 of Proposition 12. 2[ the fact that 
under the weak limit 

linnnf [Var(/i£(t, x)) — \x\e] > Var(/i(t,x)) — \x\, (3.3) 

and for x ^ 0, As{\x\^) —>■ m{m — l)|x|'"~^/2. Combine the upper bound in 
f l2Trp with identity ([23]), let £ \ in ([23]) and use Fatou's Lemma. 

To prove the lower bound of (ll.24p . let t > Cq be fixed and choose a non- 
negative smooth function f{x) with compact support such that f{x) > \x\"^ 
for \x\ < At'^^^. We have 



f{x)S{t, dx) = \im I f{x)Se{t^x)dx 
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and furthermore 



f{x)Se{t, x)dx > / \x\'^Se{t,x)dx. 
J\x\<Af2/:i 

Choose 6 > such that m + 6 < 3. By Chebyshev's inequahty and Theorem 

[ \xrSeit,x)dx < A-h-^^/^ [ \xr+^Seit,x)dx < A-^Cf"'/^ 

Hence, for appropriately chosen A, 

f{x)S,{t,x)dx > {C-^/2)e"'/\ 

Since this is true for all such /, we conclude that the lower bound of (11.241) 
holds. I 

Proof of lil.23\) . We are unable to do this by direct approximation due to 
lack of control of moments of hs{t,x) higher than 2. By direct calculation 
E[hs{t,x)] = t/4\ and we can take the e \ limit by uniform integrability 
that follows from the boundedness of VaT{hs{t, x)) argued in the beginning of 
this section. Consequently 

E[h{t,x)] =t/4\. (3.4) 

From 

Var(/i(t, x)) = E[{h{t, x) - h{0, 0) - t/Alf] 

we deduce 

As VaT{h{t, x)) = ^-2^[( MM+^)+MM) _ /i(o, 0) - t/4!) {h{t, x + S)~ h{t, x)) 
_ ^ Ht,.)+Ht,.^s) _ 0) - t/4!) {h{t, x) - hit, X - 5))]. 

Since increments are mean zero and stationary in space (part 5 of Proposition 
11.11) . the latter is equal to 

5-^E[{h{S), 5-x)- /i(0, -x)) {h{t, 6) - hit, 0))]. 

Define the "tent function" '{)six) = (6^^ — 6^^\x\)l\x\<5. We have shown that 

AsVa.T{h{t,x)) = (r,^5,5(t)) 

where the angle brackets denote integration and Tx(ps{y) = Vs{y — x) and 
we used the definition fll.211) of Sit,dx). Integrating against a test function 
i/j e C~(M) gives 

Asipix) Var(/i(t, x))dx = {ip * fs, Sit)) 
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with * denoting convolution. Let 5 \ 0. Since Var(/;,(t, x)) is locally bounded 
we can take the limit on the left. In the limit we obtain 

I j dli){x)\Bi{h{t,x))dx = {i),S{t)). 

This completes the proof of Proposition II. 61 I 

We now complete the proof of Proposition 11.61 and Theorem 11.21 with the 
following 

Proposition 3.1. For a; G M and t > 0, 

/•oo 

Var(/i(t,x)) =x + 2 / {z - x) S{t,dz). (3.5) 

J X 

This proposition implies the remaining parts of Theorem 11.21 because sym- 
metry implies 

VaT{h{t, x)) = \x\+2 {z- \x\) S{t, dz). 

J\x\ 

From this follow Yax{h{t,x)) — |x| > and the identities 



Var(/i(t,0)) = / |x| S{t,dx) 
Jr 

and for 1 < m < 3 

m(m-l) j \xr-^[Var{h{t,x)) -\x\]dx = 2 J \xrS{t,dx). 

Then we can apply the bounds from (11.240 . 
Proof of Proposition First a lemma. 

Lemma 3.2. Suppose v is a continuous, symmetric function on R and in the 
sense of distributions v" = 2/i for a symmetric probability measure /i on R. 
Assume that Jj^\x\ dfi < oo. Then 

v{x) = \x\ + v{0) - [ \z\ ^{dz) + 2 [ {z- \x\)ij{dz). (3.6) 
Jr J\x\ 

Proof. Suppose first that v" /2 is a continuous probability density. Then 

POO 

g{x) = {z — x)v"{z) dz 



satisfies g" = v" and thereby 

v{x) = ax + b+ I (z — x)v" (z) dz 
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for constants a, b. From symmetry deduce a = 1. Taking x = identifies 
b = f (0) — I J^\z\v"{z) dz. Now (13. 6p holds for smooth v. Take a symmetric 
compactly supported smooth approximate identity {(t>s}5>o, apply (13. 6p to 
(j)s * V and let 5 \ 0. I 

Continuing the proof of Proposition 13.11 apply ( 13.60 to v{x) = Var(/i(t, x)) 
to get 

Var(/i(t,x)) = |x| + ( Var(/i(t,0)) - / \z\S{t,dz) 



oo 



(3.7) 

+ 2 1 {z- \x\)S{t,dz). 



From the Appendix we get Var(/i(t, x)) — |x| ^0. Combining this with above 
gives first 

VaT{h{t,0))= \z\S{t,dz) (3.8) 



(3.9) 



and then 

Var(/i(t, x)) = |x| + 2 / (z - \x\) S{t, dz) 

J\x\ 

POO 

= x + 2 (z - x) S{t,dz). I 

J X 

4. Second class particle estimate 

In this section we prove the key estimate for the moment of a second class 
particle. The context is the asymmetric simple exclusion process (ASEP) 
jumping to the right with rate p = 1/2 and to the left with rate q = 1/2 + 
Throughout e G (0,1/4), with the real interest being the limit e \ 0. 
Probabilities associated to this process are denoted by when the process 
is stationary with Bernoulli p occupations. The macroscopic flux function is 
He{p) = —e^^'^p{l — p) and the characteristic speed Vf = H'^{p) = —5^/^(1 — 
2p). 

Let denote the probability measure of the basic coupling of two processes 
C^(^) < C(^) with this initial configuration: C~(0,0) = < 1 = C(0,0), and 
for X 7^ 0, (~{0,x) = C(0)^) have mean p and they are independent across 
the sites x. Let x(t) denote the position of the discrepancy between C~(^) and 
C{t), in other words, the position of the second class particle started at the 
origin. The mean speed of the second class particle is the characteristic speed 
(Corollary 2.5 in [2] or Theorem 2.1 in [3]): 

E?[x(t)] = Vft. (4.10) 

From (I2.10p , Theorem 12.31 is equivalent to the case p = 1/2 of the following 
theorem. 
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Theorem 4.1. With the same constants as in Theorem \l.^ for all < e < 

1/4, 1 < m < 3, and t > Coe"^, 

^^^m/3^2m/3 < |x(t) - Vftl""] < CiS'^/H^'^/^ (4.11) 

The remainder of the section proves Theorem 14. 11 with separate subsections 
for the upper and lower bound. 

4.1. Proof of the upper bound for Theorem 14.11 

Lemma 4.2. Let B G (0, oo). There exists C E (0, oo) and Ci{B) G (0, oo) 
such that the following bounds hold for all < p < 1, u > 1, < e < 1/4, and 
t>ci{B)e-^l^. 

(i) For Be^'H'^/^ <u< 20t/3, 

P^(|x(t) - Vft\ >u)< Cet\-^EP\x{t) - Vft\ + C€t\-^ + e""'/^*. (4.12) 

(ii) For u > 20t/3, 

P^(|x(t) - Vft\ >u)< e-"/^. (4.13) 

Proof. First we obtain the bounds for P^(x(t) < Vft — u). By an adjustment 
of the constant C we can assume that u is a positive integer. Fix a density 
< p < 1 and let A G (0,p). Consider a basic coupling of three ASEP's 
C > C~ ^ V with this initial configuration: 

(a) Initially {C(0,x) : x ^ 0} are i.i.d. Bernoulli(p) and C(0,0) = 1. 

(b) Initially C{0,x) = ({0,x) - 6o{x). 

(c) Initially {rj{0,x) : x ^ 0} are i.i.d. Bernoulli(A) and ^^(0,0) = 0. The 
coupling of the initial occupations is such that ({0,x) > r]{0,x) for all x ^ 0. 

Recall that basic coupling means that the processes share common Poisson 
clocks. 

Let x(t) be the position of the single second class particle between ((t) and 
C~(t), initially at the origin. Let {Xi(t) : i G Z} be the positions of the ( — rj 
second class particles, labeled so that initially 

■ ■ ■ < X_2(0) < X_i(0) < Xo(0) = < Xi(0) < X2(0) < ■ ■ ■ 

Let these second class particles preserve their labels in the dynamics and stay 
ordered. Thus the ({t) configuration consists of first class particles (the r]{t) 
process) and second class particles (the Xj(t)'s). Let P^ denote the joint 
probability distribution of these coupled processes. The marginal distribution 
of (C, C~) x) under P^ is the same as under P^. 

For X G Z, J^{t) is the net left-to-right particle current in the ( process 
across the space-time line segment from point (1/2, 0) to (x -1-1/2, t). Similarly 
J^(t) in the rj process, and J^'^it) is the net current of second class particles. 
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Current in the ( process is a sum of the first class particle current and the 
second class particle current: = J^it) + Ji'^if)- 

Basic coupling preserves x(t) G {Xj{t)}. Define the label m{t) by x(t) = 
Xm{t){t) with initial value m(0) = 0. The label m{t) performs a walk on 
the labels of the {Xj} with rates p to the left and q to the right, but jumps 
permitted only when Xj particles are adjacent. Through a comparison with a 
reversible walk, Lemma 5.2 in [3J gives the bound 

^e{m{t) < -k) < exp{-e^/^k} for alH > and A; > 0. (4.14) 

To get the first step of the estimation, note that if x(t) < Vft — u and 
m(t) > — fc, then X^k{t) < [V/tJ — u. Then among the ( — rj particles only 
X -k+i, ■ ■ ■ , Xq could have crossed from the left side of 1/2 to the right side of 
[V^^'tJ — M + 1/2 during time (0, t]. Thereby J^yP^^_^(t) < k, and by an appeal 
to (gmi) we have 

P,(x(t) < Vft -u}< P,(m(t) < -k) + P,(jf^^,^j_Jt) - J'^y.^^^^it) < k) 

< exp{-e'/'k} + P,(jf^^,^j_ Jt) - Jly.,^_Jt) < k). (4.15) 

We work on the second probability on line (14.151) . In stationary density 
p E^[Jx(t)] = Hs{p)t — px for X G Z. Process ( can be coupled with a 
stationary density p process C^^^ with at most one discrepancy. In this coupling 
I'^xi'^) ~ Jx^'^ (^)l — 1 ^'^d so we can use expectations of stationary processes at 
the expense of small errors. Let ci below be a constant that absorbs the errors 
from using means of stationary processes and from ignoring integer parts. It 
satisfies |ci| < 3. 

^ej[y.,^_St) - ^eJly.,^_St) = m{p) - mip) - U)p - tH,{X) 

+ {tH'^ip) -u)X + ci 
= -\tH:{p){p-\f + u{p-X) + c, 
= -te^l'^{p-\f + u{p-\) + ci. (4.16) 
Case 1. Be^/H'^'^ <u< bpe'^'H. Note that 20t/3 > ^pe^/H. Choose 

A = p - A.ue~^/H~^ and k = [^u^e-^/H-^\ - 3. (4.17) 
By assuming t > C{B)e^^^'^ we guarantee that u>l and 

lye-^'h-' > l^Bh'/h^/'' > 4. 

Then 

k > ^u'e-^'H-^ > 4. (4.18) 

In the next inequality below the —3 in the definition f l4.17p of k absorbs ci from 
line (I4.16p . Let X = X — EX denote a centered random variable. Continuing 



18 M. BALAZS, J. QUASTEL, AND T. SEPPALAINEN 

with the probabihty from hne (14.151) : 



< Ceeu-' (Var, [j[y.,^_M + Var, [Jly.,^_M ) ■ (4-19) 

C is a constant that can change from hne to hne but is independent of aU 
parameters. 

We develop bounds on the variances above, first for J'^ . Utihze the couphng 
with a stationary density p process. Then apply the basic identity 

Var^[j,.(t)] =p(l -p)E^|x(t) -xl (4.20) 

that links the variance of the current with the second class particle. (This is 
proved in Corollary 2.4 in [2] and in Theorem 2.1 in [3].) We find 

Vare[jf,,^P,j_(t)] < 2Var,^[jLv^P,j_„(t)] +2 

= 2p(l - p)E,|x(t) - [V/tJ +u\+2 

< E^|x(t) -IZ/tl +4u. (4.21) 

For the second variance on line (14.191) we begin in the same way: 
Var,[j;V/*j-«W] <2Var^[%Pij_.(t)] +2 

< 2A(1 - A)E^|x^(t) - [VPt\ +u\+2 

< E^|x^(t) - Vft I + 4m. (4.22) 

Here we switched to a stationary density A process and introduced a second 
class particle x'^ in this process. In order to get the same bound as on line 
(I4.2ip we wish to switch from x'''(t) to the second class particle x(t) in the 
density-p process. To this end we utilize a coupling developed in Section 3 of 
[3] . Because the density p process has higher particle density than the density 
A process, the second class particle in density A moves on average faster in the 
direction of the drift. Theorem 3.1 of [3J allows us to couple x'^ and x so that 
x(t) > x'^(t) with probability 1. Thus continuing from line f l4.22p . 

Var, [J'^y.,^_M < [x(t) - x\t)] + E,^|x(t) - Vft\ + Au (4.23) 

Now E^ [x(t) - x^(t)] = {VP - V^)t = 2e^'H{p - A) and from the choice flCTI) 
of A, 2e^/'^t{p - A) < u, hence 

Var, [jy.^^M < E^I^W - Vft\ + 5u. (4.24) 
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Insert bounds flCTD and (Km into (ICTll to get 

Pe(4v^.,j_„(t) - J^y^.i.Jt) <k)< CEt\-^W,\^{t) - + Ceeu-\ 

(4.25) 

Insert <KJ^ and fOSD into line fliTTSD to get 

P,(x(t) < V^t -u)< Cet\-^W^\i^{t) - + Ceeu-"" + 6^"'/^°* (4.26) 

and we have verified (I4.12P for P^(x(i(:) < V^t — u) for Case 1. 

Case 2. u > hpe^f^t. Let Zt be a nearest-neighbor random walk with rates 
p = 1/2 to the right and q = 1/2 + e^/^ to the left. We have the stochastic 
domination < x(i:) because no matter what the environment next to x(t), it 
has a weaker left drift than Zt. Then, since Vf = -£^/'^{l-2p), 2pe^/H < 2m/5, 
and e < 1/4, 

P^{x(t) < VPt -u}< P{Zt < -e^^H - fw}. (4.27) 
For a e (0, 1], utilizing (e° + e-°)/2 < 1 + and e~" > 1 - a, 

E,[e-"^*] = exp(-(l + e'^^)t + t^-^^{l + 2e'^^) - e'/He^'^ ) 
< exp(a2t + (a + 2a2)£i/2t). 

We can estimate P{Zt < — e^/^t — !«} < exp(— + 2a2t) and choose 
a = 1 A ^ to obtain 

Pa^it) < Vft -u}< H[-^o-'^"^ - < 20t/3 

' ^ ~ [exp(-3M/10) u>20t/3. ^ ' 

Combining (g^SD and (OSD gives Lemma for P^(x(t) < - m). 

The corresponding upper tail bound P^(x(t) — V^t > u) is obtained from 
that for P^(x(t) — V[t < —u) by a particle-hole interchange followed by a re- 
flection of the lattice. For details we refer to Lemma 5.3 in [3j. This completes 
the proof of Lemma 14. 2[ I 



Proof of the upper hound of Theorem \4-l\ Integrate Lemma|12]to get the bound 
( 14. lip on the moments of the second class particle. First for m = 1. 

POO 

E^|x(t) - Vft\ = / P^{ |x(t) - Vft\ > u} du 
Jo 

< lCB-^EP\^{t) - Vft\ +(^B+ £1/3^2/3 

Ci{B) is a new constant that depends on B. Set B = C^/^ to obtain 

E,^|x(t) - Vft\ < |Cl/3£l/3^2/3 ^ C'itl/3£-l/3exp{-J-tl/3^2/3| ^ |)g-VC_ 
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We can fix a constant cq large enougli so that, for a new constant C, 

E^|x(t) - Vft\ < Ce^^H^I^ provided t > coe~^ (4.29) 

Restrict to t that satisfy this requirement and substitute (I4.29P into Lemma 
14. 2[ Then upon using u > Be^^^t^^^ and redefining C once more, we have for 

Be^H'^/^ <u< 20t/3: 

P^( |x(t) - VPt\ >u)< Cefu-^ + 2e-"'/^*. (4.30) 
Now take 1 < m < 3 and use fl4.30p together with Lemma 14.21 

POO 

E^|x(t) - Vftr = m PP{ |x(t) - Vft\ > uju""-^ du < ^'"^Ws^sm/s 

poo poo 

+ Cmet^ / {u""-^ + 2me-"'/^*M'"-^) du + 2m e""/^?/""^ du. 

JB£l/3t2/3 i20i/3 

This gives E^|x(t) - Vftl"" < -^e'^/s^am/s provided t > co£"^ for a large 
enough cq. I 

4.2. Proof of the lower bound of Theorem 14.11 By Jensen's inequality it 
suffices to prove the lower bound for m = L Let Cub denote the constant in 
the upper bound statement that we just proved. We can also assume Cq > 1. 
Fix a constant 6 > and set 

ai = 2CuB + 1 and as = 8 + V32b + 8 a/CVb- 

Increase b if necessary so that 

b'^-2a2>l. (4.31) 

Fix a density p E (0,1) and define an auxiliary density A = p — bt^^^^e^^^^. 
Define positive integers 

u = [ait^/s^i/sj and n = [VH\ - [VH\ + u. (4.32) 

By taking cq large enough in the statement of Theorem 14. II we can ensure that 
A e (p/2,p) and m G N. 

Construct a basic coupling of three processes r] < i]^ < ( with the following 
initial state: 

(a) Initially t] has i.i.d. Bernoulli(A) occupations {t]{0,x) : x ^ n} and 
7]{0,n) = 0. 

(b) Initially r]^{0,x) = r]{0,x) + 5n{x) for all x G Z. x'^"^(t) is the location 
of the unique discrepancy between 77 (t) and r}~^{t). 

(c) Initially ( has independent occupation variables, coupled with 77(0) as 
follows: 

(c.l) ({0,x) = ri{0,x) for < X < n and C(0,'^) = 1- 

(c.2) For X > n and x < variables C(0,2:) are i.i.d. Bernoulli(p) and 

aO,x)>v{0,x). 
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Let P denote the probability measure of the coupled processes. Label the 
C — rj second class particles as {X^it) : m G Z} so that initially 

■ ■ ■ < X_i(0) < < Xo(0) = n = x(")(0) < Xi(0) < X2(0) < ■ ■ • 

Let again the random label m(t) satisfy x.^"\t) = Xm(t)it), with initial value 
m(0) = 0. In basic coupling m(-) jumps to the left with rate q and to the 
right with rate p, but only when there is an X particle adjacent to Xm{-)- As 
in the proof of the upper bound, Lemma 5.2 in [3] gives the bound 

P{m(t) >k}< exp{-e^/^k} for alH > and A; > 0. (4.33) 

By the upper bound already proved and by the choice of ai, 

P{x(")(t) < [VH\} = P{x(")(t) <n + [Vh\ - u} 

ri +2/3^-1/3 

< M-iE|x(")(t) -n- [VH\ I < , < 



CuBt'^'e'/' ^ 1 (4-34) 

[ait2/3el/3J - 2- 



This gives a lower bound for the complementary event, 

I < P{xW(t) > IVH\} < P{m(t) > A:} + P{jf^,,j(t) - Jf^p^j W > 'k}- 

(4.35) 



The reasoning behind the second inequality above is as follows: x(")(t) > [VPt\ 
and m(t) < k imply Xk(t) > [VH\ and consequently —k < J^yp^^it) = 

•^lVPt\ (^) ~ "^[VPtl (^) • 

Put k = [aat^/^e^/^J - 2. Observe from K33^ that P{m{t) > k} < e'^ < 
1/4 follows from a2t^^^e^^^ > 2e~-^/^ + 3, which is guaranteed by t > Cq^"^ and 
the definition of 02. Hence 

i < P{JlV.*j(^) - J'lvniit) > -a2t'"e'/' + 2} 

< P{JlWj(^) > -2a2t'/'£'/' - te^'\2p\ - A^) + 1} 
+ P{ J[Vptj W < -^2^'/'^'/' - te'/'(2pA - A^) - 1}. (4.36) 

Consider line fl4.36p . The rj process can be coupled with a stationary P^- 
process with at most one discrepancy. The mean current in the stationary 
process is 

E^[hvn\{t)] = tH{X) - X[VPt\ > tH{X) - XV^t = -te^'^{2pX - A^). 
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Hence 

line ( Km < P^{JivH\{t) < -a^t^'h^'^ - te^/\2p\ - A^)} 

< P^{lwn\{t) < -a^t'/'e'/'} < a^^t-^/^e-^/^ Var^ [V^^j (t)] 
E^\x{t)-[Vn\\ E^\^{t)-VH\ 2b 1 

< CuBa^' + 1^ + ^ < |. (4.37) 

After Chebyshev above we applied the basic identity (14.201) for which we intro- 
duced a second class particle x(t) in a density A system under the measure P^. 
Then we replaced \y''t\ with VH and applied the upper bound and properties 
of 02. 

Put this last bound back into line fl4.36p to get 

I < P(^) where A = {JivPt\ (t) > -2a2t^^h^^^ - te^'^{2p\ - A^) + 1} 

(4.38) 

Let 7 denote the distribution of the initial C(0) configuration described by (a)- 
(c) in the beginning of this section. As before is the density p i.i.d. Bernoulli 
measure. The Radon-Nikodym derivative is 

Bound its second moment: 

= -p{^^ ^T^)" - P"'^"^'"'^'^"^'""^ ^ -2(p). (4.39) 

Here condition t > co£:~^ implies a bound C2(p) < c>o independent of t and e. 
From f l4.38p and Schwarz's inequality 

i<P(^) = I P-{A)i{duj) 

= J < c,ipy/'{P^iA)Y^\ (4.40) 

Note the stationary mean 

E'[Jivpti{t)] = -te'/^p^ + pVH - pYVH\ < -te^'^p" + 1. 
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Continue from line f l4.40p . recalling f l4.3ip : 

(64c2(p))-^ < P'^iA) = P^iJwnm > -2a2t^'\^"' - te^'\2p\ - A^) + 1} 

= P'{Jivni{t) > {b' - 2a2)t'/'e'/'} < P'{Jivni{t) > t'/'e'/'} 
<r2/V^/='Var^[JLyptj(t)] < t-^/''e-^/^EP\^{t)-VPt\. 
This completes the proof of the lower bound and thereby the proof of Theorem 

5. Appendix: Properties of the solution 

Throughout this section h{t, x) = — log Z{t, x), where Z{t, x) is the solution 
of (11.51) starting with a two sided Brownian motion {B{x) : x G M} with 
-B(O) = 0. The goal is to prove 

'Va.T{h(t,x)) — \x\ ^ as |x| — >• oo. (5.1) 

Define the current of m = dxh across x up to time t by 

N{t,x) = h{t,x) - h{0,x) 

and the mass of u in the interval [0, x] at time t by 

M{t,x) = h{t,x) - h{t,0). 

Proposition 5.1. Vai{h{t,x)) - \x\ = Cov{N{t,0), N{t, x)) . 

Proof. Since h{0, 0) = we have h{t, x) = M{t, x) + N{t, 0). From the invari- 
ance of white noise (5 of Prop [TTT]) . Var(M(t,x)) = \x\. Hence 

Var(/i(t, x)) - \x\ = Vai{N{t, 0)) + 2 Cov(M(t, x), N{t, 0)). (5.2) 

We claim that 

Cov(M(t, x), N{t, 0)) = -Var(A^(t, 0)) + Cov(M(t, x), N{t, x)) 

+ CoY{N{t,x),N{t,0)). ^^'^^ 
To see this, note that we always have the conservation law 
N{t, x) - N{t, 0) = M(t, x) - M(0, x). 

Hence 

Cov{M{t,x),N{t,0)) = CoY{M{t,x),N{t,x)) 

- Cov(M(t, x), {M{t, x) - M(0, x))). 

But Cov{M{t,x),{M{t,x) - M{0,x))) = iVar(M(t, x) - M(0,x)). By the 
conservation law again |Var(M(t, x) — M(0,x)) = |Var(A^(t, x) — N(t,0)). 
Finally, by the translation invariance, |Var(A^()f:, x)—N(t, 0)) = Var(A^(t, 0)) — 
CoY{N{t,x),N{t,0)). This gives fOl) . 
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From fl5.3l) we can rewrite the right hand side of (15.21) as 

Cov(iV(t, x), N{t, 0)) + Cov(M(t, x), N{t, 0) + N{t, x)). (5.4) 

The proof is completed by noting that the second term vanishes by symmetry. 
To see it, note that Var(/i(t, —x)) = Var(/i(t, x)) and by translation invariance 
Cov{N{t, -x), N{t, 0)) = Cov{N{t, x), N{t, 0)). Hence Cov(M(t, x), N{t, 0) + 
N{t,x)) = Cov{M{t,-x),N{t,0) +N{t,-x)). But translating by x gives 
Cov(M(t, -x), N{t, 0) + N{t, -x)) = Cov(-M(t, x), N{t, x) + N{t, 0)). I 

Proposition 5.2. lim|^|^oo Cov(A^()f:, 0), A^()f:, x)) = 0. 

The two propositions combine to prove (15.10 . 

The proof of Proposition 15 . 21 is based on the following lemma. We need some 
notation. Fix R> and let Wi(t, x), W2(t, x) be cylindrical Wiener processes 
and Bi{x),B2{x) two sided Brownian motions with i?i(0) = -82(0) = 0, cou- 
pled as follows: For any ip G C^(]R) supported in {—oo,R), {ip,Wi(t)) = 
{(p,W2(t)) are independent of / ipdBi = f ipdB2, while for any ip G C^(M) 
supported in (_R, 00), {ip,Wi(t)), {ip,W2(t)), f ipdBi and / ipdB2 are indepen- 
dent. We will say that (Wi, dBi) and {W2, dB2) are the same on (— 00, i?) and 
independent on (i?, 00). 

Lemma 5.3. Let Zi{t,x), i = 1,2, be the solutions of U.5\) with Wi, i = 1,2 
and initial data Zi{0,x) = exp{_Bj(x)}, where (Wi,dBi) and (W2,dB2) are 
the same on {—oo,R) and independent on {R,oo). Then there is a finite C 
such that for i? > |x| -|- 2t, 

E[{Zi{t, x) - Z2{t, < Ce-^+^^'+l'l). (5.5) 

Proof. Let p{t, x) = be the heat kernel. We can write 



Zi{t,x)= pit,x-y)e'''^y^dy- / pit - s, x ~ y)Ziis,y)W,idsdy). (5.6) 




First we obtain a preliminary bound on E[Zf{t,x)]. By Schwarz's inequality 
it is bounded above by (dropping the i for clarity), 

2E[{ jpit, X - y)e''^y^dyf] + 2E[{ j j^p{t-s,x- y)Z{s, y)W{dsdy)f]. 
By Jensen's inequality, 

E[{ j p{t, x-y) e^v{B{y)]dyf] < j p{t, x - y)e^^y^dy. 
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Furthermore, 



E[{ I ^ p{t-s,.-y)Z{s,y)W{dsdy)y] 

-^^L=p{t-s,x-y)E[Z'^{s, y)]dsdy. 




Call g{t,x) = E[Zf(t,x)] and let Pt denote the heat semigroup. g{0,x) = e 



2\x\ 



and have shown that 

9{t)<Pt9{0)+ f -^P,_,g{s)ds. 

J 

Iterating once we see that this is bounded above by 



(1 + sv^)Pi(7(0) 



/ ]. , Pt— s 



Y'7r(s-u) 



Ps-ug{u)duds. 



(5.7) 



(5.8) 



The last term can be simplified by noting that Pt-gPs 
Fubini's theorem, and using 



^(t-s)(s-u) 



TT. The result is 



g{t) < (1 + s^/t/^)Ptg{0) + [ Pt-sg{s)ds. 

Jo 



Ps-u, applying 



(5.9) 



If we let g{0) = g{0) and g{t) satisfy (15. 9p with equality instead of inequality, 
then (7-51 satisfies {g-g){t)> Pt- sig - g) (s) ds with {g-g){0)=0. By the 
maximum principle for the heat equation, g < g. g{t,x) is readily computed 
with the result that for some finite C, 



S[Zf(t,x)] = (?(t,a;)<Ce^(*+l^l). 



(5.10) 



By (15. 6p again we have /(t, x) := E[{Zi(t, x) — Z2it, x))^] is bounded above 
by twice 



E 



+E 



Jpit,x- y){exp{Bi{y)} exp{B2{y)})dy 



2i 



(5.11) 
(5.12) 



Jj;^p{t-s,x~y){Z,{s,y)W,{dsdy) - Z2{s,y)W2{dsdy)) 

Explicit computation gives that (15.111) is equal to 

p{t,x- yi)p{t,x- y2)E[{e^^^y^^ - e^2fei))(eBi(y2) _ e^^fc))]^^^^^^ 




R J R 
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By Schwarz's inequality, flS.llI) < 2 J p{t,x — y)\(^y>R^e^ydy. Another explicit 
computation gives that fl5.12l) is equal to 

pt pR pt poo 

/ / p^{t-s,x-y)f{s,y)dyds + 2 / / p^{t-s,x-y)E[Z'f{s,y)]dyds 

Jo J-oo Jo Jr 

J pit-s,x-y){f{s,y) + l{y>RyCe^^-'+\y\^)dyds. 

Hence /(t) satisfies the same equation as g{t) in ( 15. 9p except that this time 
/(0,x) < l{y>R}e^^. The same argument now shows that there is a finite C 
such that f{t, x) < Ce--f^+^(*+l^l) for R>\x\+ 2t. I 

Proof of Proposition \5.2[ Let us use the notation N[t,x) for the normalized 
current N{t, x) — E[N{t, x)]. First of all note that by (11.191) and 5 of Propo- 
sition [LT], E[N^{t,x)] < C{t) and does not depend on x. Now let {Wi.dBi), 
(W2,c?-B2), iWz,dB-i) be coupled so that iWi.dBi) and (W^2, 0^-82) are the 
same on (— oo,a;/2) and independent on (a;/2,oo), (W2,dB2) and (W3,dB3) 
are the same on {x/2,oo) and independent on {—oo,x/2), and (Wi,dBi) and 
(W^ydBs) are independent. Let Ni, N2, N3 be the currents corresponding to 
the three different pairs. Of course Cov(iV(it, 0), N{t, x)) = E[Ni{t, 0)Ni{t, x)]. 
By Schwarz's inequality 

\E[Ni{t, 0)iVi(t, x)] - E[N2{t, 0)iVi(t, x)] \ < C\E[{N,{t, 0) - iV2(t, 0))Y^^ 

(5.13) 

and 

\E[N^{t,0)N2{t,x)] - E[N^{t,0)N3{t,x)]\ < C\E[{N2{t,x) - N,{t,x))^^^^^^ 

(5.14) 

By independence, E[Ni{t,0)N3(t,x)] = 0. By symmetry, 
E[N2{t, 0)N,{t, x)] = E[N,{t, 0)N2{t, x)]. 

Hence 

Cov(iV(t, 0), N{t, x)) < C{E[{h,{t, 0) - /i2(t, 0))2])i/2. (5.15) 

For each L > 0, let log^, 2; = logz for z > and log^, 2; = — logL for 
< z < We have 

E[(/ii(t,0)-/i2(t,0))2] < 2£^[(log^Zi(t,0)-logiZ2(t,0))2] 

+4E[l|o<Zi<L-i}(logZi(t,0))2]. (5.16) 

Because log^^ is Lipschitz with constant L we have 

E[{\ogL Z^it, 0) - log^ Z2(t, 0))2] < L2E[(Zi(t, 0) - Z^it, 0))2]. (5.17) 

By Lemma for each fixed L, this vanishes as |x| 00. On the other hand, 
since £'[(log Zi(t, 0))^] < 00, by the dominated convergence theorem, 

lim E[l|o<z,<L-i}(logZi(t,0))2] = 0. (5.18) 

L — ^■oo 




SCALING OF KPZ/STOCHASTIC BURGERS 

This completes the proof. I 



27 



References 

J. Baik, P. Deift, and K. Johansson. On the distribution of the length of the longest 
increasing subsequence of random permutations. J. Amer. Math. Soc, 12:1119- 1178, 
1999. 

M. Balazs and T. Seppalainen. Exact connections between current fluctuations and the 
second class particle in a class of deposition models. J. Stat. Phys., 127(2), 2007. 
M. Balazs and T. Seppalainen. Fluctuation bounds for the asymmetric simple exclusion 
process. ALEA Lat. Am. J. Probab. Math. Stat, 6:1-24, 2009. 

A.-L. Barabasi and H. E. Stanley. Fractal concepts in surface growth. Cambridge Uni- 
versity Press, Cambridge, 1995. 

Lorenzo Bertini and Giambattista Giacomin. Stochastic Burgers and KPZ equations 

from particle systems. Comm. Math. Phys., 183(3):571-607, 1997. 

S. Bczcrra, S. Tindcl. and F. Vicns. Supcrdiffusivity for a brownian polymer in a 

continuous gaussian environment. Ann. Probab., 36(5):1642-1675, 2008. 

P. Billingsley. Convergence of probability measures. Wiley, 1968. 

Terence Chan. Scaling limits of Wick ordered KPZ equation. Comm. Math. Phys., 

209(3):671-690, 2000. 

P. L. Ferrari and H. Spohn. Scaling limit for the space-time covariance of the stationary 

totally asymmetric simple exclusion process. Comm. Math. Phys., 265(l):l-44, 2006. 
Dieter Forster, David R. Nelson, and Michael J. Stephen. Large-distance and long-time 
properties of a randomly stirred fluid. Phys. Rev. A (3), 16(2):732-749, 1977. 
H. Holden, B. 0kscndal. J. Ub0c, and T. Zhang. Stochastic partial differential equa- 
tions. A modeling, white noise functional approach. Birkhauser Boston, Boston, 1996. 
K. Johansson. Transversal fluctuations for increasing subsequences on the plane. 
Probab. Theory Related Fields, 116:445-456, 2000. 

K. Kardar, G. Parisi, and Y.Z. Zhang. Dynamic scaling of growing interfaces. Phys. 
Rev. Lett, 56:889-892, 1986. 

H. Krug and H. Spohn. Kinetic roughenning of growing surfaces, pages 412-525. Cam- 
bridge Univ. Press., 1991. 

C. Licea, C. Newman, and M. Piza. Supcrdiffusivity in first-passage percolation. Prob. 

Th.Rel Fields, 106:559 591, 1996. 

O. Mejane. Upper bound of a volume exponent for directed polymers in a random 
environment. Ann. Inst H. Poincar'e Probab. Statist, 40:299-308, 2004. 
C. Mueller. On the support of solutions to the heat equation with noise. Stochastics 
Stochastics Rep., 37(4):225-245, 1991. 

David Nualart and Boris Rozovskii. Weighted stochastic Sobolev spaces and bilinear 
SPDEs driven by space-time white noise. J. Funct Anal, 149(1):200 225, 1997. 
M. Piza. Directed polymers in a random environment: Some results on fluctuations. J. 
Statist Phys., 89:581-603, 1997. 

M. Prahofcr and H. Spohn. Current fluctuations for the totally asymmetric simple 
exclusion process. Progress in Probability. Birkhauser, 2002. 

G. Da Prato and J. Zabczyk. Stochastic equations in infinite dimensions. Cambridge 

University Press, Cambridge, 1992. 

J. Quastel and B. Valko. t^^^ supcrdiffusivity of finite-range asymmetric exclusion 
processes on F. Comm. Math. Phys., 273(2):379-394, 2007. 



28 



M. BALAZS, J. QUASTEL, AND T. SEPPALAINEN 



[23] J. Walsh. An introduction to stochastic partial differential equations, volume 1180 of 
Lecture Notes in Mathematics, pages 265-439. Springer- Verlag, 1986. 

Department of Stochastics, Budapest University of Technology and Eco- 
nomics 

E-mail address: balazs@math.bme.hu 
URL: www.math.bme.hu/~balazs 

Departments of Mathematics and Statistics, University of Toronto 
E-mail address: quastel@math.toronto.edu 
URL: www .math . toronto . edu/quastel 

Department of Mathematics, University of Wisconsin-Madison 
E-mail address: seppalai@math.wisc.edu 
URL: www. math. wise. edu/^seppalai 



